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Abstract: We show that in the absence of a magnetic field the spectrum of the 
magnetic Schrodinger operator in an annulus depends on the cosine of the flux 
associated with the magnetic potential. This result follows from an analysis of a 
singularity in the "wave trace" for this Schrodinger operator, and hence shows that 
even in the absence of a magnetic field the magnetic potential can change the asymp- 
totics of the Schrodinger spectrum, i.e. the Aharonov-Bohm effect takes place. We 



also study the Aharonov-Bohm effect for the magnetic Schrodinger operator on a 



torus. 

51. Introduction. 
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Let O be the exterior of a bounded region in M 2 with smooth boundary, and let 

£3 ■ I i 

H A>V = -(-id Xl + A^x)) 2 + ~(-id X2 + A 2 {x)) 2 - V(x) (1) 



This is the Schrodinger operator for a particle of mass 1 and charge -1 moving 
in O under the influence of the magnetic potential A = (Ai, A2) and the electric 



C^*- ■ potential V. We assume that 

2: d X2 A x - d Xl A 2 = in n, (2) 

i.e., the magnetic field vanishes in O. Given a simple, closed curve 7 in O encircling 
the complement of O, we define the magnetic flux by 



^ 1 Q! 7 = / A(x) ■ dx. 

5-H 



In view of (2) ct 7 only depends on the orientation of 7. 

In the seminal paper [AB] Aharonov and Bohm showed that if a 1 7^ mod 27r, 
then one can detect the cosine of the magnetic flux in the scattering of particles 
in this quantum system, i.e. the magnetic potential has a physical impact even 
when the magnetic field is zero in O. This effect is called the Aharonov-Bohm 
effect. They computed the scattering cross-section explicitly for O = M 2 \{0}, when 
A(x) = (— X2/\x\ 2 , xi/\x\ 2 ) and V(x) = 0. They also proposed an experiment to 
demonstrate this effect. However, the first generally accepted experimental verifi- 
cation of the Aharonov-Bohm (AB) effect was done many years later by Tonomura 
et al. [T] . For further mathematical work on the AB effect see [N] , [W] , [RY] , [E2] , 
[E3], [EIO]. 
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In [H] Helffer showed that A(x) can influence the spectrum of Ha,v when the 
magnetic field is zero in 0. In the semi-classical setting with V(x) — > oo, as \x\ — > oo, 
and O = {\x\ > 1} he showed that the lowest Dirichlet eigenvalue depended on 
the cosine of the magnetic flux. Earlier related results on magnetic Schrodinger 
operators are due to Lavine and O'Carroll ([L-C]). 

In this paper we study the Schrodinger operator in the domain = Q D {\x\ < 
R} with Dirichlet boundary conditions on \x\ = R and dVt. We compute the 
singularity at t = 3R\^3 of the distribution trace of the fundamental solution of the 
initial-boundary value problem 

uu + Ha,vu = 0, u(x, 0) = f(x), ut(x, 0) = 0, u(x, t) = when x G 8Qr. (3) 

This distribution trace is known as the "wave trace" for this problem, and it is 
given by 

CO 

^cos(t\/A~), 

i=i 

where {Xj}j? =1 are the Dirichlet eigenvalues of Ha,v in Hence its singularities 
are determined by the behavior of the Xj as j — > oo. These singularities are well- 
known to appear only at the lengths of periodic broken ray paths in Qr. The 
singularity at t = 3R\^3 comes from equilateral triangles in with vertices on 
\x\ = R. To compute this singularity we need to know that 3R\/3 is isolated in 
the set of lengths of broken periodic rays. To ensure that we assume that the 
complement of O, fi c , is strictly convex and contained in {\x\ < 1} and R>8 (see 
Remark 1.1), but any assumption that makes the length of the inscribed equilateral 
triangles isolated in the lengths of periodic reflected ray paths will suffice. The 
geometry that we have chosen makes the singularity unchanged when one changes 
the sign of ct 7 . Hence we cannot recover more than the cosine of a 1 from it (see 
Remark 1.2). 

A definitive computation of leading singularities in wave traces was given by 
Duistermaat and Guillemin in [DG] for manifolds without boundary. For manifolds 
with boundary the analogous computation has not been done in that generality. To 
carry it out in here we have taken this opportunity to present a different method of 
computation that replaces Fourier integral operators with superpositions of Gauss- 
ian beams (cf. [CRR] and Chapter 5 of [CR]). The result of the computation is 
that the leading singularity at t = L = 3Ry/3 is the real part of the distribution 

±^R{2R)VW cos(^ A(x) ■dx) (t _ L 1 _ . Q)3/2 , 

where the initial ± does not depend on A. 

In the final section of this paper we consider Ha,v on (flat) 2-torus and obtain 
essentially the same result: under a nondegeneracy assumption on the torus the 
singularities in the wave trace at times equal to the lengths of curves in a homology 
basis determine the cosines of magnetic fluxes around those curves (see Theorem 
5.1). 
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Remark LI: The only fact from geometry needed here is (and we only need it 
for circles): rays and their reflections inside an ellipse are tangent to the same 
ellipse confocal with the boundary ellipse. So rays in \x\ < R tangent to a circle 
|x| = r > 1 will never enter \x\ < 1 after reflection in \x\ = R, while rays that enter 
| a; | < 1 will always re-enter |x| < 1 after reflection in |x| = R. Since the boundary 
curve C is convex, rays entering \x\ < 1 will leave \x\ < 1 after at most one reflection. 
This gives the following bounds on the length L of periodic ray paths that hit C. 
For rays that close after entering \x\ < 1 k times 

2kR-2k < L < 2kR + 2k. 

So periodic rays that enter \x\ < 1 more than three times have lengths are greater 
than 8R — 8, and the equilateral triangles are the (isolated) shortest periodic rays 
that never enter \x\ < 1 (assuming R > 2). So we need AR + 4 < 3R\^3 < 6R — 6. 
That happens as soon as R > 8 (picking the first whole number that works). 

Remark 1.2: If O = > 1} and V = 0, the mapping u(x) — > u(—x) sends 
eigenfunctions of Ha,o to eigenfunctions of H-a,o bijectively. Thus the wave traces 
of these operators must be identical. The leading singularity in the wave trace 
at t = 3V3R does not depend on the boundary of O or V(x), hence it will be 
unchanged when A is replaced by —A in these cases, too. Therefore, one cannot 
distinguish ct 7 and — a 7 using the leading singularity. The same ambiguity arises 
in the results in [AB] and [H]. 

§2. Singularities of the Wave Trace. 

Let E(x,y,t) denote the fundamental solution for the initial-boundary value 
problem (3) . The wave front set of the distribution kernel of E is contained in the 
canonical relation for the bicharacteristic flow (see Melrose-Sjostrand, [MS I, II]. 
For this problem this canonical relation is defined as follows: Let u(x) denote the 
outer unit normal to 8Qr at x. Given (yo,r}o) with yo e and \rjo\ = 1, define 
(x(s,y,r)),£(s,y,r])) = (y + sr],r)) until, at s = si, yi = x(s 1 ,rio,y ) E dQ, R . Then, 
if rj ■ v{yi) ^ 0, continue (x(s, y , 770), £(s, yo, Vo)) for s > s i as (yi + sr)i,Vi), where 
Vi = Vo ~ 2(i/(yi) ■ rjo)iy(yi). Continue the bicharacteristic this way, reflecting when 
x(s, yo, rjo) hits 90^, as long as x(s, yo, rjo) does not intersect dVt^ tangentially. At 
points of tangential intersection one has to distinguish grazing and gliding points. 
However, since we assume that the boundary of O c is strictly convex, points of 
tangential intersection with dVt are grazing points and bicharacteristics continue 
unaffected by these intersections. When y is in the interior of Qr, a bicharacteristic 
with initial data (yo,Vo) wih never intersect \x\ = R tangentially. Hence, the wave 
front set of the kernel of E(-, -,t) is the union over yo G and r]o G S 1 of the 
points 

(x(t, y , Vo),£,(t, ?/o, %), ?/o, -%), 

where (x(t, yo, rjo), £(t, yo, rjo)) are the reflected bicharateristics described above. 
Strictly speaking, the wave front set is the closure of that set and includes a "bound- 
ary wave front set" over \x\ = R (see [MS] for details). 
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Since E(x,y,t) is a distribution in t depending smoothly on (x, y) e x 
J n E(x,x,t)dx is well-defined, and we have the following relation 



The singular support of T is contained in the set of t such that (yo, ?7o, 2/o, —Vo) £ 
VF-F(i?(x, y, £)) for some yo G H_r, [GM]. The choice of Q and here implies that, 
for t in a sufficiently small neighborhood of 3Ry/3 : (yo, ?7o, yo, — Vo) £ WF(.E(a;, y, £)) 
only if the ray x(s, yo, ?7o) traces an inscribed equilateral triangle. 

To compute the singularities in the wave trace we need a parametrix for the 
initial-boundary value problem (3). Since this parametrix will differ from E(x, y, t) 
by an integral operator with a smooth kernel, we can use it to compute sin- 
gularities. Since we are only interested in singularities arising from inscribed 
equilateral triangles, we only need a parametrix which captures the singularities 
of f E(x,y,t)f(y)dy when WF(f) C {y, rj) : y E fl R ,\y ■ rj^] = R/2}, where 
(771,772) = (772 , — ?7i ) - These singularities hit 8Qr nontangentially, and hence this 
parametrix construction can be done with reflection at the boundary. This obser- 
vation applies equally well to constructions with Fourier integral operators and the 
Gaussian beam superpositions used here. 

§3. The Gaussian beam construction. Here we will outline the construction of 
a parametrix for (3) , for initial data with wave fronts projecting onto the inscribed 
equilateral triangles. We will continue to let rj have length one. The Gaussian beam 
method allows one to do the following (see [R] for more details): 

i) For any ray, (x(t),t) = (z + tr/,t), in space-time, one can construct a function 
4>(x,t; z,r]) satisfying: 

(a) For any given integer N, ((fit) 2 — |<^x| 2 vanishes to order N on (x(t),t) and 
Iva{(fi xx } is positive definite on (x(t),t). 

(b) <fi(x, 0; z, rf) = x ■ i] + ^\x — z\ 2 on \x — z\ < 5, and (fit(x, 0; z, rj) = —1. 

Moreover, if V is a curve with unit normal v at x(to) and rj is not tangent to T, 
then one can construct (fi r = (fi on T, satisfying (a) for the reflected ray (x(t ) + (t — 
to)rj r , t), where rf = u — 2(u ■ rj)u. 

ii) Once (fi has been constructed, for any given integer N, one can solve the transport 
equations 



2(fi t (a J ) t -2(fi x -(a J ) x + 2iA(x)-(fi x + (fi tt -A(fi)a j = -(d? - (d x +iA(x)) 2 )a J _ 1 , j >0 



to order iV on (x(t),t), and impose the initial conditions ao(0, x; z, rj) = 1 and 
aj((0, x; z, rj) = for j > on \x — z\ < 5. 




2<fit(a )t - 2(fi x ■ (a )x + (2iA(x) ■ (fi x + <fi u - A</>)a = 0, 



(4) 
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For the singularity computation we need to know the leading amplitude ao on 
the ray beginning at z in direction 77. 

We define a(x, t; z, 77, r) to be the formal sum 

a(x,t;z,r],r) = ^a,j{x,t; z,rj)r~ J . (5) 

j>0 

As before one can reflect in a plane curve V which is transverse to the ray, and we 
impose a r = — a on T to satisfy Dirichlet boundary conditions. 

Using the preceding constructions we can construct the operator 
[V(t)f](x)= 1 -([V + (t)f](x) + [V_(t)f](x)), 

where 
[V±(t)f](x) = 

V— i- f e ir ^ {x > ±t > z ^ a k (x,±t;z,ri,r)f(rri)r 2 drdr,dz. (6) 

( 27r ) JR + xS 1 x{\z\<R+S} 

Here, 0° is the phase function with 0; 2, 77) = x • rj + ^\x — z\ 2 , and for k > 0, 

e ir ^ ( *> t; ^V(2;,t;£,?7,r) 

is the (Dirichlet) reflection of e ir ^ fc 1 ( x,t ' z,v ^a k ~ 1 (x,t; z,rj,r) in the circle |x| = -R. 
Since Gaussian beams can be constructed to for any finite ray segment, we can 
assume that each term in (6) is defined on < 2R} when necessary. Note that 
in this notation the variables (z, if) in (j) k remain the initial data at t=0 for the 
ray where Im{</> fc } = 0. Note also that the integration in r in (6) is in the sense of 
distributions. 

For the parametrix construction we need V^(0)/ = f+Kf where K is an operator 
with a smooth kernel. From (6) we have 

[V(0)f]( x ) = — ^ / e irx -^ x - y \ 2/2 f(rri)r 2 drdridz. 

(27TJ Jr+xS 1 x{\z\<2R\ 



Since 



(27T)* 



f e irx -^ x - y ^ /2 f(rri)r 2 drdridz = f(x) 

JR+xS 1 xR 2 . 



and / is supported in {\x\ < R}, it follows that omitting the contribution from 
{\z\ > R + 5} in (6) only adds an operator with a smooth kernel. 

To compute singularities of the wave trace we need to make the kernels of the 
operators V±(t) explicit. The distribution kernels of these operators are sums of 
terms of the form 



S(t)= [ e ir(t> ^ z ^- irv - y a{x,t;z,ri,r)r 2 drdridz, 

JR + xS 1 xR 2 z 



(7) 
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As was stated earlier, these operators are smooth in (x,y), and we can compute their 
traces by integrating these kernels over the diagonal y = x. Thus the (distribution) 
trace of V(t) is a sum of terms of the form 

Tr(cf>,a)= [ e^^^-^aix^t^z^^ydrdrjdzdx. (8) 

We want to compute the singularity in t of this trace at t = L = 3R\/3, and we 
only need to consider t in \t — L\ < 5, where 5 is small enough that {t : \t — L\ < 6} 
contains no other lengths of periodic rays in the disk \x\ < R. 

§4. Calculation of the singularity at t = L = 3\^3R 

For rj = (?7i , 772 ) with \rj\ = 1 define rj = (772, — ^i), the "right hand" normal. 
To compute the singularity at t = L we only need the parametrix restricted to 
R/2 — e < \z • rj | < R/2 + e for any fixed positive e. Since the broken ray 
x(t, z, rj) is initially of the form x = z + trj, rj • z > corresponds to rays going 
counterclockwise around z = 0, and rj 1 - ■ z < corresponds to rays going clockwise 
around z = 0. 

In the preceding section we concluded that the singularity in the wave trace at 
t = L could be calculated from a sum of integrals of the form 

1^2 f Q r2 dr f dV (/ a o(^ 5 ±i ^^) e "^ (x ' ±t ' Z ' ?7) " X ' ?7) ^^}- ( 9 ) 

The integral in r is to be taken in distribution sense. Until the end of this section 
we will consider (9) in the case that the phase <p is the beam phase resulting from 
reflecting the bicharacteristic with initial data (x, £) = (z,rj) three times in \x\ = 
R. The amplitudes ao(x,t, z,rj) are determined by the transport equation (4). 
The contributions to the singularity from the + and — terms in (9) are complex 
conjugates of each other, and from here one we only consider the "+" term. 

We assume that that ao vanishes when \z ■ rj \ is not close to R/2. Note that 
we can assume that (f)(x, t, z, rj) is defined for all (x, z, t) when | - 77 j is sufficiently 
close to R/2. 

The main step in isolating the singularity is an application of the method of 
stationary phase to (9). For that we introduce the change of coordinates 

x = u + vrj + wrj 1 - , z = vrj + wrj 1 - , u G IR 2 , v, w G K. 

Our objective is the elimination of the integral in (w, w) by stationary phase. To 
see when the phase is real and stationary in these variables note that 

i) the phase is real only when x = x{t, z, rj), 

ii) the derivative of the phase with respect to u at x = x(t, z, rj) is 

<Px ~rj = £{t,z,rj) - rj, 

which vanishes precisely when three reflections have made £ return to its initial 
value. That implies \z • rj | = R/2. Since the reflected ray will return to z when 
t = L and it is propagating in the direction rj, x{t,z,rj) = z + (t — L)rj. Hence 

11. = (t. — T,\n and \m\ = R/9. nn thp. sta.tinna.rv sfit in ?/, 
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The derivative of the phase with respect to w at x = x(t, z, rj) is 

v ± ■ <Px + v ± ■ <Pz - v ■ v ± 

which vanishes, since (f) z (x(t, z, rj), t, z, rj) = (f) z (x(0, z,rj),0, z,rj) = d z (x ■ rj + i\x — 
z\ 2 /2)\ x=z = 0. Thus we will need to do the stationary phase computation at 
(u,w) = ((t-L)r,,±R/2). 

Calculation of asymptotics by stationary phase requires the computation of the 
determinant of the Hessian of the phase, and here this computation is rather long. 
We have found it useful to consider the phase and the bicharacteristics defined for 
all rj 7^ by homogeneity. That makes the Jacobian matrix 



--de} 



a b 
c d 



symplectic. Using <f> x (x(t, z,rj),t, z,rj) = £(t,z,rj) and (f) z (x(t, z,rj),t, z,rj) = and 
setting M = (f) X x(x(t, z, rj), t, z, rj), one computes directly that at x = x(t, z, rj) 



H = 



def 



^zx 



'xz 
■'zz 



M 



Ma 



c l - a l M a l Ma - a l c 



Letting O v be the matrix with columns rj and rj , one sees that the Hessian of the 
phase in (9) with respect to the variables (it, v, w) is B t HB where 



B=\ I °- 

' O v 



However, we need the Hessian with respect to (u,w). We will see that 



is a 



null vector for H, and we have B 




1 

Vo/ 



Moreover, letting P v denote the 



orthogonal projection of M 2 onto (rj) , one computes 



"S P. 



B = 



/0 0\ 



10 

Vo o o o/ 



Hence, 



det 








'u 2 u 1 YU2U2 u 'ru 2 w 
10 

fiwui ^wuz (ftww 



UlUI ^ 

/ 



, M c-Ma \ , . 

(U1 1 c'-a'M a'Ma-a'c + Pj - [W) 
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To proceed with this computation we need to know F(t). The computation 
begins with the formulas for x(t, z, rj) and £(£, z, rj) after three reflections: 

X(t, Z, rj) = Wj^ + (t + - Q\/R 2 - w 2 ) ||| 

and, setting rj = |??|(cos6>,sin6>), 

IX) IX) 

£(t,z,rj) = |t7|(cos(0 + 7T - 6sin _1 — ), sin(0 + n - 6 sin -1 — )). 

R R 



One checks that d z w = ^ and d v w = — (z ■ an d this implies that the 

Jacobian §| at w = ±R/2 is ^\rj\P v ±. So c = ^|r?|P^. Using d v 6 = -r?V|r/1 2 , 
one finds that at w = ±R/2 

dv~ v v R V V R \v\ v 

So d = I- ^vP v ±. 

The computations of the derivatives of x(t, z, rj) are longer, but they are simpli- 
fied by the observation that z,r))\ = \rj\. At w = ±R/2 one has 



So a = /+(£- L + v)^P v ± 

dx 



To compute j^- at w = ±R/2 one uses 



at to = ±-R/2, and the less obvious result that 

/£\ / 1 . 4=V3z-rj rj rj 



Combining those with d v v = (z ■ 77 )^3 = rk-^^rj , one has 



_i_ 
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dx Tj 
drj 



\m 



■(-(z-rj) 



T] 



_L 



\m 

Thus 

F(t) = 



+ o((± =F 2y/3(z -v)^,-) + (t-L+ ^)(] 

z l'/l ' 



4^/3 z -J] 



\m 



(t - L + t3) 4^3 



(1- 



z • r\. 



R 77 



R \r]\ 



)P^)\v\~ 1 



(z ■ rj) 



2 ^V 1 - • 



I+(t-L + v)^P v ± 

R \'l\- r V ± 



(t-L) 



1 - 



4V3 
R 



v )p . _ 4V3 «ip 



(11) 



Now we can resume the computation of the Hessian. First we compute the deter- 
minant of the Hessian. For this the only facts that we need from the computation 
of the symplectic matrix F(t) - it is a good check on the computation to verify that 
it is symplectic - are that a, 6, c and d commute with P v with aP v = dP v = P v 
and bPrj = bP v = 0. We will also eventually use the exact form of c. Note that 
since F(t) is symplectic a l c and d l b are symmetric and a l d — c l b = I. 

Returning to (10) we have 



and 



M c-Ma \ ( I a\ = ( M c 

c l - a l M a 1 Ma - a*c + P v J \0 I J ~ ^ c* - a*M P v 

I 0\ f M c\ _ f M c \ 

Since M = (c + id)(a + ib)' 1 (cf. [CRR]), 

M c \ ( a + ib 0\ _ ( c + id c 
c* a*c + P v J \ / J ~ \ c f a + ic% etc + P v 

-a* I\ ( c + id c \ _ /i(c*&-a*d) P v \ _ ( -il P v 

I / V c*a + ic l b etc + P„ J \ c + id c ) \ c + id c 



Finally 



-ic + d I\ ( -il P v \ _ ( P v +c 
I J \c + id c J ~ \-U P v 



From the preceding, using the exact form of c, one can read off the determinant of 
the Hessian of the phase (at u = (t — L)rj, w = ±R/2). It is 

(-l)(^)det((a + z6)" 1 ). (12) 
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We do not need to know the Hessian matrix for the phase, but we will need to 
know the null space of its imaginary part. Note that H ^ = ^o) ( no ^ e ^ na ^ 

a f c = c*a). This says that the columns of ^ ^ are in the null space of lm{H}. 

Since Im{M} is positive definite, the null space of Im{H} can be at most two 

dimensional. Since the columns of (^j^j are linearly independent, the nullspace 

of Im{H} must be the span of those columns. Hence from (11) one sees that the 
nullspace of lm{H} has the orthogonal basis 

{(^((l + tt-L + ^^hW)}. (13) 

At this point it is convenient to calculate the amplitude ao- Note that 4>t{ao)t — 
4>x ■ (flo)x = — ~dt a o( x (ti z, v)i t-> z i v) = def — a-Q. Hence (4) implies that, after three 
reflections, evaluating ao when z = vrj ± (R/2)rj ± gives 

a (x(t,Z, V ),t,Z, V ) = ( _ 1) 3 e ^>(x( S ))x( S )d Se (/ o t [^-A^(x( S ), S ) ds )/2_ (M) 

Note that \(f) x \ + 4> t vanishes to second order on (x, t, z, rj) = (x(t, z, rj), t, z, r])), and 
on that ray (j) t t + 4>tx • x = and <p x = rj/\r}\. Hence, differentiating \4> x \ + 4>t = 
with respect to x and using <$> tt = —4>t x • rj/\r}\, we have 

(ptt -A(j> = r}/\r}\ ■ Mr}/\r}\ - trace(M). 

Differentiating x = with respect to z and rj one sees that a + ib = (I — P v )(c + 
id). Hence, using M = (c + id)(a + z6) _1 

log det(a + ib)) = trace((d + ib)(a + ib)' 1 ) = trace(M - P V M) 

= trace(M - P^MP V ) = A<f> - (ptt- 
Thus, we can conclude that after three reflections 

a (x(t,z, V ),t,z, V ) = (-l) 3 (det(a + z6))- 1 /V/o M*(*))*(s)ds ^ (15) 

where the square root is chosen continuously in t starting at 1 when t = 0. We have 
A(x(s))x(s)ds = a 7 , where 7 is the equilateral triangle traced by x(s, z, rf) with 
z = vrj + (R/2)r] 1 - or z = v rj — ( J R/2)?y ± . Since the magnetic field vanishes in O, a~ 
is independent of v and rj, and its value when z = vrj + (R/2)rj 1 - is the negative of 
its value when z = vrj — (R/2)r) ± . 

Now we can evaluate the integral in (u, w) asymptotically by the method of 
stationary phase. From here on we will always have rj = (cos 9, sin 9). The standard 
form of the stationary phase lemma, ([Hor], Theorem 7.7.5), gives the following: if 
f(y) is a smooth function such that Im{/} > 0, f y (yo) = and the Hessian f yy (yo) 
is nonsingular, then for a smooth with support in a sufficiently small neighborhood 
of 2/0 7 one has the asymptotic expansion 



/ 



e irf ( y) a(y)dy = — ^ c,/- 

\ T / j=Q 
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and the leading coefficient is given by 

Co = e ir ^a(y )(det(-if yy (y ))- 1/2 . (16) 

The subscript * indicates that the square root has to be chosen according to a rule 
determined by the eigenvalues of f yy , see [Hor]. 

In our case we will use stationary phase to eliminate the integrations in u and 
w in (9) - recall that z = vrj + wrj and x = u + vrj + wrj . The stationary point 
Uq in (13) is either (it, w) = ((t — L)rj, R/2) or (it, w) = ((t — L)rj, —R/2) . Since 

<p(x(t, z, rj), t, z, rj) = <p(x(0, z, rj), 0, z, rj) = z ■ 77, 

and we have 

f(Vo) = <P(x(t, z, rj),t, z, 77) - x(t, z,7])-r] 

evaluated at (it, if) = ((£ — L)rj,R/2) or (u,w) = ((t — L)rj, —R/2), it follows that 
f(yo) = — (t — L). The domain of integration in (it, v, w, rj) is 



{(it, v, w, rj) : \rj\ = 1, |it + vrj + wrj 1 - \ < R and \/ w 2 + v 2 < R + 5}. (17) 

We consider (9) as an iterated integral with the integrations in (it, w) done first. 
After we use the stationary phase lemma in those integrations, the resulting inte- 
grand is evaluated at (u,w) = ((t — L)r),±R/2), and, since we can assume that 
\t — L\ is smaller than 5, the domain of integration in (y, rj) becomes 

D = def l-^R -(t-L),^R-(t- L)] x S 1 . 

There is a problem when v is at the endpoints of its range: then the restriction 
I it + vrj + wrj 1 - 1 < R does not allow integration in (it, w) over a neighborhood of the 
stationary points, and the stationary phase lemma does not apply. We will handle 
this by first extending the domain of integration in (ix, w) to a neighborhood of those 
boundary points, and then showing that the integration over the added portion of 
the domain does not contribute to the leading term in the singularity. Assuming 
that the domain of integration has been enlarged, (12), (15) and (16) give uniformly 
for (y, rj) E D 

[ a (x,t,z,rj)e ir ^ x ' t ' z ^- x ^ ) dudw = ±^^K(t)e~ irit - L) +0(^), (18) 
Jd(v, v ) r 6/2 r b ' 2 

where D(v,rj) = {(u, w) : \u + vrj + wrj \ < R + e} for e small, and c(R) = 
(27r) 3 / 2 (^j) 1 / 2 e~ ;t7r/ ' 4 . The choice of sign ± is determined by (15) and (16): it is 
-1 when the square roots of det(a + i6) implicit in (15) and (16) agree and +1 when 
they do not. The factor 

K(t) = exp(if A(x + (s)) ■ x + (s)ds) + exp(z / A(x~(s)) • x~(s)ds) 
Jo Jo 
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arises from adding the contributions from stationary points with w = —R/2 and 
w = R/2. The path x~ (s) with w = —R/2 goes clockwise around the origin, and the 
path x + (s) with w = R/2 is counterclockwise. Hence K(L) = 2cos(/ 7 A(x) ■ dx), 
and integration over (y, r\, r) gives the leading singularity in the trace at t = L 
announced in earlier, i.e. 

±n 3 R(2R) 1 / 2 3 1/4 cos( f A(x) ■ dx)(t — L — zO)" 3 / 2 . (19) 

This uses the distribution calculation 

-*<*-*)r r i/2 dr _ e-^T(3/2) 
(t-L-iO)W 



L 



The final step in this argument is showing that (19) really is the leading term 
in the singularity. Hence we need to look at the integral integral (in (it, v, w)) over 
the region added to derive (18). This region is given by 

fo /a 

E = def {(u, v,w): R< lu+vrj+wr] 1 - \ < R+e, v E [-^—R-(t-L), -^—R-(t-L)}}. 

Note that, since the stationary points are given by (it, w) = (t — L)rj, ±i?/2), the 
stationary points are in E only when v is at one of the endpoints of its range. 
Without loss of generality we assume that v is the endpoint ^R— (t — L). We will 
only work on a small neighborhood of that point, since the domain only needed to 
be enlarged when v was in a neighborhood of that point. Likewise we will ignore 
all contributions beyond the second order Taylor expansion of the phase, and use 
the tangent plane as the boundary, since higher order corrections contribute lower 
order terms in r to the integrals. 

To estimate the integral in (it, w) we introduce an orthogonal frame, 
{ei, e2, 63} with e\ an inner (with respect to E) normal to |x| = R at 

(x, Z) = (x , Zq) = (u + V 7] + Wo^^VoT] + W^) 

A convenient choice for this frame is 

ei = >/3(t7, rj) + (V 1 , V 1 ), e 2 = Vs(v, ~v) + (v^ ■> and 

e3 = (?7, -v) ~ v 7 ^?^, ~V ± )- 

Note that e 2 and es are orthogonal to (77, rj) which spans the null space of the 
Hessian of the phase. Using this frame we can introduce coordinates on R 4 with 
origin at (xq, zq) by setting 



(x -xq,z- z ) = u\e\ + u 2 e 2 + (v - v )(r], rj) + we 3 . 
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Since e 2 and e 3 are orthogonal to (rj, rj), it follows from (12) that the Hessian of 
the phase at (xq, zq) with respect to u 2 and w is nonsingular. Since those variables 
are tangential to the boundary |x| = R at (xq, zq)- remember that we are replacing 
boundaries by tangent planes - the integral with respect to those variables is order 
r _1 by stationary phase. Moreover, after this use of stationary phase the phase 
function is evaluated at (x, z) = (xo, zq) + u\e\ + (v — vo)(rj, v)- Hence, (see (13)) 
the imaginary part of the phase modulo terms which are O(-uf) is given by au\ 
with a > 0. Now we are left with integration in (v, ui). The domain of integration 
in these variables is determined by the constraint \x\ > R on E. Letting s = v$ — v 
(note that s > 0) this constraint translates to 

and it follows that for (ui, s) sufficiently small we must have u± > as for some a > 
(any a < V3/4 will work). Thus the leading integral in (r, £ti, s) has the form 

pOO /»CXD pOO 

I(t)= / / re- ir ^ t - L ^ ai+ia ^ d2 ^ 1 ^(p(u 1 ,v)drdu 1 ds, 

JO J as JO 

where d 2 > Cu 2 , <p(ui,v) has a compact support (p(ui, v) = 1 near (0, 0), ct2 > 0. 
Integrating in r we get 

/•OO />oo -j^ 

I(t) = / 7 — 7 -, : — \ 1 „ n ip(ui,v)duids. 

Jo Jas [-(t-L) + (a 1 + i a 2 )d^ K ' 

Note that \{t - L) + (a t +ia 2 )d 2 \ 2 > d(\t - L\ + uj) 2 > C t \t - L\ 1+e \u\ for 
< e < 1/2. Then 2 - 2e > 1 and 



J as 



du\ 1 



u\- 2 " (l-2e)(asy-^- 



' as "i 

Therefore estimating the integrals in u\ and s we get 

Thus I(t) has a lower order singularity near t = L. 

One should check that the contribution from the reflected beam near xq is also 
of lower order. Recall that the reflected phase (jf described at the beginning of §3 
equals <fi on the boundary. Hence its Hessian with respect to the variables (u2, w) at 
(xq, zq) is identical to the corresponding Hessian of <f>. Moreover, the reflection rule 
implies that the tangential component of cjf x is equal to the tangential component of 
4> x on the boundary. This implies that the phase 4> r — rj ■ x is stationary at boundary 
points when — rj • x is stationary, and we can use the stationary phase argument 
in the preceding paragraph again to eliminate the integrals in (v,2,w). 
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This leaves an integrand for the integration in ui which is order r -1 as before. 
However, since u\ is normal to the boundary the resulting phase is not stationary 
in u\ and integration by parts in u\ makes the contribution from the reflected beam 
of lower order. This completes the proof that (19) is the leading singularity in the 
wave trace at t = L. 

§5. The Aharonov-Bohm Effect on a Torus. 

The Aharonov-Bohm effect only arises when the underlying domain is not simply 
connected. In the previous sections the domain was an annulus. Here we consider 
the Schrodinger operator on a torus. Let L = {miei + m 2 e 2 '■ m G Z 2 }, where 
{ei,e2} is a basis for M 2 . We assume that the lattice L has the property: For 
d, d' G L, if \d'\ = then d' = ±d. This is a generic condition that implies that 
the group of isometries of L consists of lattice translations and the inversion d — > —d. 
Associated to L one has the dual lattice L* = {5 G K 2 : 5 ■ d G Z for all d G L}. 

We consider the Schrodinger operator from (1), 

H A ,v = \{~id Xl + A 1 (x)) 2 + \(-id X2 + A 2 (x)) 2 - V(x), 

acting on functions on T 2 = R 2 /L. The functions A = (Ai, A 2 ) and V are assumed 
to be smooth on T 2 and hence they have smooth extensions to M, 2 satisfying A{x + 
d) = A(x) and V(x + d) = V(x) for all d G L. As before we assume that the 
magnetic field vanishes 

a i2 A 1 -a xi A 2 = 0onT 2 . (20) 
Thus for any closed curve 7 on T 2 the flux 

ct 7 = / A(x) ■ dx, 

is determined by the homology class of 7. We let 71 and 72 be a basis for the 
homology group, for instance 

7j = {te J ,tG[0,l)}, J = 1,2, (21) 

and denote the corresponding fluxes by a\ and a 2 - 

Let g{x) G C°°(T 2 ) be such that \g(x)\ = 1. The conjugation of Ha,v by the 
unitary operator of multiplication by g(x) transforms H^ y to v , where A = 
A+ig~ 1 Vg. The condition \g(x) \ = 1 on T 2 implies that g(x) = exp(2ni5-x+ip(x)), 
where 5 G L* and <p(x) is periodic. Hence a\(A) = a\(A) — 2tvS ■ e\, a 2 (A) = 
012(A) — 2tvS ■ e 2 . Therefore if A and A are gauge equivalent we have 



OLj{A) = otj(A) modulo 2-7T, j = 1,2. 



(22) 
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Expanding A(x) in a Fourier series we have 

A(x) = A + J2 A 5^ l5 ' x , 

6eL*\{0} 

where A = |T 2 | _1 J J2 A(x)dx, |T 2 | denotes the area of {sei + tei\ < s, t < 1}. 
Since d X2 Ai = d Xl Ai we have A(x) = Aq + V(p(x), where 

seL*\{0} 

Therefore when (20) holds A(x) is gauge equivalent to the constant potential Aq. 
Two constant magnetic potentials Aq and Aq are not gauge equivalent if (22) does 
not hold. When Aq is not gauge equivalent to either Aq or — Aq the potentials Aq 
and Aq have a different physical impact, in particular, the spectra of Ha ,v and 
H i T/ are not the same. 

The last assertion is a consequence of the following theorem. 

Theorem 5.1. Suppose (20) holds. The spectrum of EL Ay as a self-adjoint oper- 
ator on L 2 (T 2 ) determines coscti and cos«2, where ctj = A(x) ■ dx, j = 1,2. 

Theorem 5.1 complements the results of [G], [ER1] and [El]. In particular it 
shows that, if A and A give rise to zero magnetic fields on T 2 but different values 
for coscti and coso^, the Schrodinger operators, Ha,v and v will have different 
spectra. This proves the Aharonov-Bohm effect on the torus. 

Proof of Theorem 5.1. As in the preceding sections we start with the wave trace 
formula 

00 „ 

cos(t\/~\j) = / Ej2(x,x,t)dx, 

where {A-,}"^ is the spectrum of Ha,v on T 2 and Ej2(x,y,t) is the solution to 
E tt + H A yE = on T 2 x R satisfying E(x, y, 0) = 5(x - y) and E t (x, y, 0) = 0. 
Note that 

E T 2(x,y,t) = ^2e r2 (x + d,y,t), 
deL 

where E R 2 is the solution to E tt + Ha,vE = on R 2 x R satisfying E(x, y, 0) = 
8{x — y) and E t (x,y,0) = when Ha,v has been extended to R 2 by making its 
coefficients periodic, i.e. A(x + d) = A(x) and V(x + d) = V(x) for all deL. Hence 

/ Ej2 (x, x, t)dx = y / E-g2{x + d, x,t)dx. (23) 
M deL Jji 

Since E^ 2 is smooth off the cone \x — y\ 2 = t 2 , and our assumption on L implies 
that only two lattice vectors can have \d\ 2 = t 2 for a fixed value of t, the singularity 
in the wave trace at t = \d\, must come from (cf. [ERT], [ER2]) 



/ E R2 (x + d,x,t)dx+ / E-g 
Jt 2 Jt 2 



! (x — d, x, t)dx. 
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To compute the leading singularities in this trace we will use the Hadamard- 
Hormander parametrix (cf. [Hor]). We have 

E R 2(x,y,t) = d t (E+(x,y,t) - E+(x,y, -t)), 

where E + is the forward fundamental solution. The Hadamard-Hormander para- 
metrix construction for E + writes E + as an asymptotic sum of terms with increasing 
regularity. The first term is cio(x,y)eo(\x — y\,t), where 

eo = -^-j=(t 2 — \x — y\ 2 )j^^ 2 when t > and eo = when t < 0, and 

a (x, y) = exp(z / (x - y) ■ A(y + s(x - y))ds). 
Jo 

Therefore (cf. [ER1]) the singularity of the trace at t = \d\ determines 1(d) +I(—d) 
where 

1(d) = / exp(z / d ■ A(x + sd)ds)dx. 



IT 2 JO 

Since A(x) = Aq + V(p(x), where <p(x) is periodic, we have 

/ d ■ A(x + sd)ds = d • Aq since / d • V(p(x + sd)ds = 
Jo Jo 

Therefore 1(d) = e ld ' A ° |T 2 | and hence the singularity of the wave trace at t = \d\ 
determines cos(Ao • d) for all d G L. In particular, when d = ej and 7^ = {tej, t G 
[0, 1)},J = 1,2, we get ctj = f A(x) ■ dx = ej ■ Aq. Thus the singularities of the 
wave trace when t = \ej\ determine cosctj for j = 1,2. When V(x) = V(—x), 
then Ha () ,v and H-a ,v are isospectral and one can only recover cosaj, j = 1,2, 
from the spectrum. When V is not even, the question of whether one could recover 
exp(ictj), j = 1, 2, from the spectrum is open. 
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